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THE  STRUCTURE  OF  THE  NEAR  TIP  FIELD  DURING 


TRANSIENT  ELASTODYNAMIC  CRACK  GROWTH 
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Abstract 

The  process  of  dynamic  crack  growth  in  a  nominally  elastic  material  under  conditions  of  plane 
strain  or  plane  stress  is  considered.  Of  particular  concern  is  the  influence  of  the  transient  nature  of 
the  process  on  the  stress  field  in  the  immediate  vicinity  of  the  crack  tip  during  nonsteady  growth. 
Asymptotically,  the  crack  tip  stress  field  is  square  root  singular  at  the  crack  tip,  with  the  angular 
variation  of  the  singular  field  depending  weakly  on  the  instantaneous  crack  tip  speed  and  with  the 
instantaneous  stress  intensity  factor  being  a  scalar  multiplier  of  the  singular  field.  However,  for  a 
material  particle  at  a  small  distance  from  the  moving  crack,  the  local  stress  field  depends  not  only 
on  instantaneous  values  of  crack  speed  and  stress  intensity  factor,  but  also  on  the  past  history  of 
these  time-dependent  quantities.  A  representation  of  the  crack  tip  field  is  obtained  in  the  form  of 
an  expansion  about  the  crack  tip  in  powers  of  radial  coordinate,  with  the  coefficients  depending  on 
the  time  rates  of  change  of  crack  tip  speed  and  stress  intensity  factor.  This  representation  is  used 
to  interpret  some  experimental  observations,  with  the  conclusion  that  the  higher  order  expansion 
provides  an  accurate  description  of  crack  tip  fields  under  fairly  severe  transient  conditions.  In 
addition,  some  estimates  are  made  of  the  practical  limits  of  using  a  stress  intensity  factor  field 
alone  to  characterize  the  local  fields. 


Introduction 


One  of  the  cornerstone  concepts  of  linear  elastic  fracture  mechanics  is  the  role  of  the  stress  intensity 
factor  as  a  crack  tip  field  characterizing  parameter.  Irwin  (1957)  observed  that  the  elastic  stress 
field  near  the  tip  of  a  crack  which  opens  symmetrically  (mode  I  opening)  has  the  now  familiar  form 

ff.j  =  -jj=='Zij(0)  +  0(1)  as  r  -►  0 
V27t  r 

in  a  polar  coordinate  system  centered  at  the  crack  tip.  The  dependence  on  radial  coordinate  r 
is  explicit,  and  the  angular  variation  functions  S,j  can  be  expressed  in  terms  of  trigonometric 
functions  of  angle  8.  The  leading  term  in  an  asymptotic  expansion  of  the  crack  tip  field  in  powers 
of  r  is  thus  universal.  The  only  feature  which  varies  with  loading  conditions  and  geometrical 
configuration  of  the  cracked  body  is  the  scalar  multiplier,  the  elastic  stress  intensity  factor  Kj.  For 
strictly  two-dimensional  deformations,  the  above  field  can  be  valid  only  within  a  region  near  the 
crack  tip  of  extent  small  compared  to  the  crack  length,  distance  to  the  nearest  boundary,  or  any 
other  characteristic  dimension  of  the  body.  On  the  other  hand,  the  assumption  of  a  two  dimensional 
elastic  field  cannot  be  valid  right  up  to  the  crack  tip,  even  under  plane  strain  conditions.  Within 
some  small  region  near  the  tip,  the  two  dimensional  elastic  field  must  give  way  to  a  region  of  inelastic 
deformation.  Nonetheless,  if  this  inelastic  zone  is  completely  surrounded  by  a  stress  intensity  factor 
field,  then  it  is  commonly  assumed  that  the  stress  intensity  factor  can  be  used  to  characterize  the 
fracture  process  and  that  fracture  will  begin  when  the  value  of  stress  intensity  factor  has  been 
increased  to  a  material  specific  value,  the  fracture  toughness  of  the  material. 

When  the  stress  intensity  factor  idea  is  applied  to  the  case  of  a  through-the-thickness  crack  in 
a  plate,  the  restrictions  become  somewhat  more  severe  if  the  process  is  to  be  characterized  in  terms 
of  a  stress  intensity  factor  associated  with  a  two  dimensional  state  of  plane  stress.  In  addition  to 
the  limitation  noted  above  concerning  overall  body  dimensions  and  crack  length,  such  points  must 
also  be  far  enough  from  the  crack  edge  to  be  outside  the  zone  of  influence  of  three  dimensional 
effects  around  the  crack  edge.  Typically,  the  plane  stress  assumption  is  valid  only  for  r  >  h/2 
where  h  is  the  plate  thickness  (Rosakis  and  Ravi-Chandar  (1986),  Yang  and  Freund  (1986)). 

If  the  crack  growth  process  in  a  plate  is  dynamic,  so  that  the  inertia  of  the  material  comes  into 
play,  there  is  yet  another  factor  which  complicates  the  application  of  the  stress  intensity  factor  idea 
under  conditions  of  plane  stress.  This  factor  is  due  to  the  wave  character  of  the  mechanical  fields 
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in  the  body  during  crack  growth.  To  illustrate  this  effect,  consider  a  crack  growing  dynamically  at 
a  constant  speed.  Suppose  that  at  equally  spaced  intervals  of  time,  wave  signals  are  emitted  from 
the  crack  tip.  A  set  of  such  waves  is  illustrated  schematically  in  Figure  1  where  the  circles  represent 
wavefronts  of  the  signals.  The  centers  of  the  circles  are  located  at  equally  spaced  points  along  the 
crack  growth  direction.  Consider  a  particular  field  point  P  that  is  near  to  the  crack  tip  in  some 
sense.  It  is  then  clear  from  the  diagram  that  this  point  is  far  behind  the  wavefronts  emitted  earliest, 
compared  to  the  distance  from  P  to  the  crack  tip.  On  the  other  hand,  the  point  is  outside  of  the 
wavefronts  emitted  most  recently.  If  information  measured  at  this  point,  as  well  as  at  similar  points 
around  the  crack  edge,  is  used  to  infer  something  about  the  value  of  the  stress  intensity  factor, 
then  the  wave  nature  of  the  fields  should  be  considered.  In  some  sense,  local  measurements  sense 
the  recent  history  of  the  stress  intensity  factor.  If  the  stress  intensity  factor  is  essentially  constant 
during  this  time,  then  the  effect  is  obviously  negligible.  However,  if  the  stress  intensity  factor  is 
varying  rapidly  in  time,  then  its  value  may  change  significantly  during  the  time  of  influence  at  P 
and  the  existence  of  a  region  of  practically  significant  size  around  the  crack  edge  in  which  the  field  is 
characterized  by  the  instantaneous  value  of  stress  intensity  factor  is  uncertain.  The  possibility  that 
the  extent  of  the  stress  intensity  factor  during  dynamic  crack  growth  under  transient  conditions  is 
more  limited  than  a  steady  state  analysis  would  indicate  was  first  suggested  by  the  analysis  of  Ma 
and  Freund  ( 1986). 

The  purpose  of  the  present  study  is  two-fold.  First,  a  higher  order  asymptotic  expansion  of 
crack  tip  fields  under  plane  stress  conditions  is  obtained  for  transient  crack  growth.  In  the  present 
context,  transient  crack  growth  is  understood  to  include  processes  in  which  both  the  crack  tip 
speed  and  the  dynamic  stress  intensity  factor  are  differentiable  functions  of  time.  The  analytical 
result  is  used  to  interpret  some  experimental  observations,  with  the  conclusion  that  the  higher  order 
expansion  provides  a  good  representation  of  crack  tip  fields  under  fairly  severe  transient  conditions. 
Then,  some  estimates  are  made  of  the  practical  limits  of  using  a  stress  intensity  factor  field  alone 
to  characterize  the  local  fields. 
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Governing  Equations  for  Plane  Stress  Elastodynamic  Crack  Growth 

Consider  a  homogeneous  and  isotropic  linear  elastic  solid  occupying  a  two-dimensional  region  in  the 
x'  Xjj-plane.  The  outer  boundary  is  subjected  to  traction  and/or  displacement  boundary  conditions 
of  a  type  to  ensure  uniqueness  of  solutions.  Suppose  that  a  planar  crack  grows  through  the  body 
with  nonuniform  crack  speed  u(t).  Within  the  framework  of  the  theory  of  plane  stress,  the  two- 
dimensional  displacement  vector  u  is  governed  by  the  equation  (Freund,  1990) 

e?V(V-u)-c2VA  VAu  =  ii,  (1) 

where  V  is  the  two-dimensional  gradient  operator  and  a  superposed  dot  indicates  differentiation 
with  respect  to  time.  In  terms  of  Young’s  modulus  E  and  Poisson’s  ratio  u,  the  longitudinal  and 
shear  wave  speeds  for  plane  stress  are 

c,  =  [E/ {l  -  v2)p\X/2  and  ca  =  [E/2{1  +  v)p}1/2  , 

respectively,  where  p  is  the  material  mass  density. 

Any  displacement  u  which  is  derived  from  longitudinal  and  shear  wave  potentials  <f>(x\,x'2,t) 
and  according  to 

u  =  V0  +  VA^  (2) 

with 

cfV2(f>-4>  =  0,  c2,V2^-V>  =  0  (3) 

satisfies  (1).  Conversely,  any  solution  of  (1)  has  the  representation  (2).  In  plane  stress  xp  has  a 
single  nonzero  component  denoted  here  by  ip. 

Suppose  that  a  translating  coordinate  system  ( X\,X2 )  is  introduced  with  its  origin  at  the 
moving  crack  tip  and  oriented  with  the  zpaxis  aligned  with  the  direction  of  crack  growth.  The 
displacement  potentials  can  then  be  viewed  as  functions  of  position  in  the  moving  coordinate  system, 
as  well  as  time.  The  new  functions  are  denoted  by  $  and  that  is,  $(xi,Z2,t)  =  4>{xi,x'2 ,<)  and 
'P(xi,X2,t)  =  xp(x\,x'2,t).  Under  the  transformation  of  coordinates,  the  partied  differential  equation 
(3)  governing  the  longitudinal  displacement  potential  becomes 

#2$>  /  u2\  d$  v  <9$  v  d2$  1  M  _  , .. 

dx\  \  cj  )  dx\  cj  dxi  +  cj  dxi dt  c ]  dt 2 
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The  shear  wave  potential  '®(ari,*2,t)  satisfies  the  same  equation  with  C|  replaced  by  c4.  It  is  also 
noted  for  future  reference  that  the  first  scalar  invariant  of  stress  for  the  assumed  state  of  plane 
stress  is  determined  by  the  displacement  potential  $  according  to 


<Tll  +  <*22  =  2 p(c]  -  c^)Vx2$ 


where  V,?  is  the  Laplacian  operator  in  the  translating  coordinate  system. 


(5) 


Asymptotic  Crack  Tip  Field  for  Transient  Crack  Growth 

To  derive  an  asymptotic  expansion  for  the  stress  components  as  r  =  \/x\  +  x\  — ►  0,  a  standard 
device  is  employed  whereby  the  region  around  the  crack  tip  is  expanded  so  it  fills  the  entire  field 
of  observation  (Freund,  1990).  To  this  end,  rescaled  coordinates  r)a  =  xa/e  are  introduced,  with 
e  >  0  being  a  small  parameter  and  a  =  1,2.  As  e  — ♦  0,  all  points  in  the  xi,X2-plane  except  those 
near  the  crack  tip  are  mapped  beyond  of  the  range  of  observation  in  the  77, ,  772-plane. 

It  is  assumed  that  $  has  an  expansion  in  powers  of  f  of  the  form 

OO 

$(x1,x2,<)  =  $(o?i, €772,0  =  t^$m(771, 772,0,  (6) 

m=0 

where  4>0  represents  the  main  contribution  to  the  asymptotic  solution,  $1  represents  the  first  order 
correction,  and  so  on.  The  first  term  of  this  series  (ro  =  0)  corresponds  to  the  expected  square  root 
singular  contribution  proportional  to  r-1^2  in  the  asymptotic  near  tip  stress  field. 

The  assumed  asymptotic  form  (6)  is  substituted  into  the  governing  equation  (4)  and  the 
coefficient  of  each  power  of  e  is  set  equal  to  zero  to  obtain  a  system  of  coupled  differential  equations 
for  $m.  The  first  several  such  equations  are 


202*o  d2*0 

°l  dr,\  +  dr,] 

2<92$1  32$i 

dr,*  +  O772 

2 

a‘  drj\  dr,] 

2d2$3  d2$3 

a'  dr,\  +  dr, 2 


=  0, 

=  0, 


2»1/20  f  i/2fl*o  1 
-  ■  c]  dt  1  077,  J  ’ 

_  _2 f  1/2 

c]  dt  Y  dm  J  ' 


(7a) 

(76) 


(7c) 


(7  d) 
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2^4  dH*  __2v^d_  j  1/2^j1 , 1^0 
a<  drj{  dj]\  cj  dt(  drji  J  cj  dt 2 

2d2$5  d2$s  2V1/2  d  f  1/2d$3 1  ,  1  <92$i 
° 1  drj[  drfo  c }  dt  drfi  J  cf  dt 2 

and  so  on,  where  a/  =  (1  —  v2 /c2)1^2. 

For  any  value  of  m,  the  governing  equation  has  the  general  form 


(7e) 

(7/) 


d2*m 

dr)2  dt?2 


d2$m  _  2vx  12  d  j  1/2d$m_ 2  ]  1  d2$m-4 

c?  dt  (  J  c?  dt 2 


(8) 


for  m  =  0,1,2,...,  where  =  0  if  A:  <  0.  The  shear  wave  potential  $  is  also  assumed  to  be 
expressible  in  terms  of  a  similar  expansion  for  e,  with  coefficients  satisfying  an  equation  similar 
to  (8)  with  ci  replaced  by  c,. 

It  is  noted  that,  for  the  special  case  of  steady  state  crack  growth,  v  =  0  and  d$m/dt  =  0  for 
m  =  0,1,2,....  In  such  a  case,  the  equations  (7)  are  not  coupled  and  each  reduces  to  Laplace’s 
equation  in  the  coordinates  771, 0^772 .  The  solution  for  this  case  is  discussed  by  Dally  (1987)  who 
attributes  the  original  results  to  unpublished  work  of  G.  R.  Irwin.  The  corresponding  functions 
<I>m  are  independent  of  time  in  the  moving  coordinate  system.  In  the  general  transient  case,  on  the 
other  hand,  the  only  uncoupled  equations  are  those  for  m  =  0  and  m  =  1.  Indeed,  as  will  be  seen, 
$0  and  have  the  same  spatial  structure  in  both  the  transient  and  the  steady  state  cases.  This 
is  not  so,  however,  for  if  m  >  1. 

An  expansion  for  the  first  stress  invariant  in  the  scaled  coordinate  system  tji,tj2  can  be  obtained 
from  (5)  and  (6)  as 


OO 

<7U  +  <722  =  2 p(cf  -  c]) 

m= 0 


'  d*f?m  d$m 
.  dnl  drft  , 


In  view  of  (8),  this  expression  can  be  rewritten  as 


<J\\  +  <722 

OO 

-  V 

’ v 2  d2$m 

2v*/2  d 

(1/2  d$m-2  ] 

1  ,  1  d2$m_4‘ 

2 p(cj  -  c2) 

- 

m= 0 

_c2  drj\ 

cf  dt 

l  dvi  J 

[  +  c2  dt2 

This  relation  does  not  involve  derivatives  of  with  respect  to  t)2. 


(9) 


(10) 
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The  Case  of  Transient  Crack  Growth  at  Constant  Velocity 


In  this  section,  a  special  case  of  transient  crack  growth  is  considered,  namely,  crack  propagation 
with  constant  crack  tip  speed  v  but  with  a  time-dependent  stress  intensity  factor.  Although  the 
crack  tip  speed  is  constant  in  this  case,  the  potential  functions  $  and  $  are  still  time-dependent, 
in  general,  and  the  right  side  of  (8)  is  not  identically  zero  as  it  is  for  steady  state  crack  growth. 
Crack  growth  phenomena  having  this  behavior  have  been  observed  by  Ravi-Chandar  and  Knauss 
*  (1987)  and  Rosakis  et  al.  (1983,1984),  among  others.  A  mathematical  model  having  this  feature 
is  provided  by  the  plane  stress  problem  of  a  crack  suddenly  beginning  to  grow  symmetrically  at 
constant  rate  from  zero  initial  length  in  a  uniform  body  oubjected  to  a  remote  tensile  field.  The 
solution  of  this  problem  was  first  obtained  by  Broberg  (1960).  As  is  evident  from  the  analytic 
solution  of  the  Broberg  problem,  which  is  considered  in  detail  in  Freund  (1990),  the  stress  intensity 
factor  is  a  function  of  time.  In  fact,  based  on  dimensional  considerations,  the  stress  intensity  factor 
varies  in  proportion  to  time  t  raised  to  the  power  one-half.  The  Broberg  solution  will  be  examined 
further  in  a  subsequent  section. 

The  equations  (7)  may  be  solved  sequentially,  and  the  first  steps  for  m  =  0  and  m  =  1  have 
been  presented  by  Freund  (1990).  For  these  cases,  the  solutions  for  symmetric  modes  of  crack 
growth  are 

=  Am{t)r”2  cos  |  (m  +  3)0; ,  m  =  0, 1  (11) 


where  rj  =  (r/f  +  ajr)l)  and  0;  =  tan-1(a(?72/7?i).  Note  that  $m(n,0j,f)  =  $m(Th,T}2,t)  by 

definition.  For  constant  velocity  crack  growth,  (rh  0()  are  translating  polar  coordinates  distorted 
from  conventional  polar  coordinates  by  an  amount  determined  by  the  crack  speed  v.  Because  v  is 
constant,  the  coordinates  (r(,0()  are  time-independent.  Time  enters  the  solution  (11)  only  through 
the  coefficients  Am(t).  By  definition,  the  coefficient  A0(f)  is  related  to  the  stress  intensity  factor 
(Freund,  1990)  by 


4jj(l)  1+aj 
°(  '  2nV2 *  A(») 

where  k"}(t)  is  the  time-dependent  dynamic  stress  intensity  factor  and  A(u)  =  4a/a,  -  (1  +  a2)2. 

The  functions  $2  and  $3  are  obtained  by  solving  (8)  for  m  =  2  and  3.  For  these  values  of  m, 
(8)  reduces  to 


2d2<Em  d$m  _  2v  d2$m- 2 

1  dr}\  drj*  cf  dtdrji 


(12) 
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The  result  of  substituting  (11)  into  the  right  side  of  (12)  and  converting  to  the  distorted  polar 
coordinates  (r/,0;)  is 

r-r2ji  m  ,  1  d$m  ,  1  ,-,1  r  t  i  (m  — 1)/2  ,\/i  /io\ 

~  ^  =  D  COS  2  (m  "  W  ’  (13) 

where  Dl  [.4*]  denotes  the  operator 


for  A:  =  0, 1 , 2, 3, . . which  reduces  to 

afcf  at 

for  constant  velocity  crack  growth.  In  the  development  to  follow,  this  operator  is  understood  to 
yield  =  0  for  k  <  0. 

Relation  (13)  holds  only  for  m  =  0, 1,2,3.  The  general  solution  for  these  values  of  m  and 
symmetric  mode  I  crack  growth  is 

<hm  =  r^~  |^lm(<)cos  i(m  +  3)0;  +  —  2^+2^  cos  2(m  “  l)ff‘  j  *  (14) 

which  also  applies  only  for  m  =  0, 1,2,3.  Note  that  (11)  is  a  special  case  of  (14)  which  corresponds 
to  m  =  0, 1.  For  m  =  4,5  the  procedure  is  continued,  in  which  case  (8)  becomes 


2d2<Em  d<t>m  _  2V  d2$m- 2  1  d2$m-4 

a‘  dr]\  dr}]  c]  dtdr] i  c 2  dt 2 


If  the  potential  is  again  viewed  as  a  function  of  the  seeded  polar  coordinates  (r/,0;),  and 
if  (14)  is  substituted  into  the  right  side  of  (15),  then  $(r;,0(,t)  must  satisfy  the  partial  differential 
equation 


Vi2$m(ri,0,,O  =  |  ^1[^m-2(0]  +  g(^l~i4)(/)]  +  ^m-4  COS  Hm  “  M 

,  D2[Am-4{t)}  i  \  (m  — 1)/2 

+  2(m-l)  C°S  S(m  -  5)9‘)  T‘ 


where 


D2[Ak}  =  Dl  [Dl[Ak]]  ,  Ak  = 


1  d2Ak 
a2c2  dt2 
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for  k  =  0, 1, —  Note  that  (16)  holds  only  for  integers  m  =  The  general  solution  of  (16) 

for  mode  I  crack  growth  is 


*m(n,0[,t)  =  r(,n+3,/2(Am(0cos  j(m  +  3)0, 

+  2(m‘+  t)  [d'M»-»(0]  +  +  ^-4]  cos  1  (m  -1)0,  (17) 

- 4(0]  „>„i 

+  “8(^-l)2  “M"*-2)9') 

for  m  =  1, . . .  ,5.  It  is  clear  that  (11)  and  (14)  are  special  cases  of  (17)  for  m  =  0, 1  and  m  —  0, 1,2,3, 
respectively. 

A  six  term  expansion  for  the  stress  invariant  as  a  function  of  position  around  the  crack  tip  can 
be  obtained  by  substituting  (17)  into  (10).  If  the  terms  in  the  series  are  written  out  explicitly,  then 
this  expansion  is 


0\\  +  cr22  3d2  ,  _1/2  ,  „  ,  2d2  , 

2p(c? =  4 cfAor‘  C08^‘  +  T  1 


{  l^~Ai  +  D' [A°]]  C°S  2^  +  Mo]  cos  t^}  r'/ 

{  %A3  +  (1_  icf) cos0'}r? 


+  { [% tAx  +  0  ”  £?) £>1[/l21  +  K1  ”  £?) 1)21/101  +  0 "  £? ) Ao] cos  *9t  (18) 

+  j£SDliA*\  +  ^  (1 "  ^2)  1,2  Mo]  +  l&fA°  COS  2^  +  9^2Z?2[Ao]cosf0' }r?/2 


+  { .^rA5  + 11 "  £f) Z)1[/l3] + ^2[Al] + (l  ~  if) Ai. cos2*' 

+  [^M*]  +  |  (l  -  £2)  ^Ml]  +  }  rl  + 

where  .CAfAfc]  was  introduced  in  (13)  and  D2[Ak]  was  introduced  in  (16).  Note  that  the  spatial 
variation  of  the  terms  in  this  expansion  with  coefficients  A *  are  identical  to  those  obtained  from  an 


asymptotic  analysis  under  the  assumption  .of  steady  state  crack  growth.  Of  course,  in  the  general 
transient  case,  the  quantities  represented  by  A*  are  functions  of  time.  Also,  Z)1[Ajt],  £)2[A*],  and 
.4*  are  normalized  derivatives  with  respect  to  time  of  the  coefficients  A*.  For  example,  for  constant 
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velocity  v  but  time  dependent  stress  intensity  factor  Kj{t),  the  quantity  Z)1^]  can  be  reduced  to 
the  first  time  derivative  of  the  stress  intensity  factor, 


Dl[A0] 


4t>(l-K>2)  dKj 
Haftfy/ 27rA(u)  dt 


Comparison  with  Constant  Velocity  Experiments 

In  this  section,  some  experimental  results  on  rapid  crack  propagation  are  interpreted  on  the  basis 
of  the  foregoing  asymptotic  analysis.  The  main  purpose  is  to  determine  if  the  crack  tip  fields 
actually  observed  in  the  vicinity  of  dynamically  growing  fractures  can  be  represented  by  the  fields 
anticipated  in  the  higher  order  asymptotic  analysis.  The  comparison  is  restricted  to  the  case  in 
which  transient  crack  growth  is  indeed  observed  to  occur  with  constant  crack  speed,  following 
fracture  initiation. 

A  full  field  optical  technique,  called  the  Coherent  Gradient  Sensing  (CGS)  method,  has  been 
developed  by  Tippur  et  al.  (1989)  for  measuring  crack  tip  deformation  fields.  This  method  can 
be  used  in  either  transmission  mode  for  transparent  materials  or  in  reflection  mode  for  opaque 
materials.  Like  the  well-known  method  of  caustics,  the  CGS  method  is  based  on  the  detection  of 
gradients  in  the  optical  path  caused  by  deformation  of  a  nominally  two  dimensional  sample.  In 
transmission  mode,  these  gradients  arise  from  changes  in  refractive  index  of  the  material  due  to 
nonuniform  strain  in  the  sample,  while  in  reflection  mode  they  arise  from  changes  in  slope  of  the 
sample  surface  due  to  nonuniform  out-of-plane  strain.  In  either  case,  if  the  deformation  can  be 
viewed  as  being  consistent  with  the  theory  of  plane  stress  in  the  region  in  which  measurements 
are  made,  then  the  experimental  information  is  related  to  the  in-plane  gradient  of  the  first  stress 
invariant,  which  has  components  d{au  +  cr22)/dx0  r°r  a  =  1.2. 

If  the  terms  in  (18)  are  differentiated  with  respect  to  x\  and  then  rearranged  in  order  to 
isolate  the  term  that  would  represent  the  stress  intensity  factor  if  the  field  were  indeed  a  square 
root  singular  stress  intensity  factor  field,  it  is  found  that 


V^Trr^2  d((Tn  +  <722  ) 
F(v)cos|0(  dxi 


h'i(t)  +  r, 


/MO 


cos 

cos  I#; 


+  /MO 


COS  1 0/1 

COS  |0/  J 


+  /M0r/ 


3/2 


COS  5  ffl 


+  rl 


+  /MO  +  /MO 


COS  \d\ 

cos  1 9 1 


+  /M0  r)n 


COS  0J 

cos  1 6 1 


(20) 


where  F(v)  —  (1  -f  a2)(a2  —  a2)/A(u)  and  Pi,..  .,Pi  are  functions  of  time  and  crack  tip  velocity 
determined  by  the  various  coefficients  in  (18).  As  a  matter  of  convenience,  the  left  side  of  (20) 
will  be  denoted  by  Yf(ri,8[,t)  in  subsequent  discussion.  The  important  observation  is  that  if,  at 
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any  particular  instant  of  time,  the  near  tip  stress  field  is  accurately  characterized  by  a  square  root 
singular  stress  intensity  factor  field  over  some  region,  then  Yf(ri,9i,t)  is  independent  of  r;  and  9i 
throu0hout  that  region.  Furthermore,  the  value  of  Yf  at  such  a  time  is  precisely  the  value  of  the 
corresponding  dynamic  stress  intensity  factor  Kj(t).  On  the  other  hand,  if  there  is  no  region  in 
the  vicinity  of  the  crack  tip  throughout  which  Yf  is  uniform,  then  it  must  be  concluded  that  the 
local  crack  ti^  field  is  not  accurately  represented  by  a  stress  intensity  factor  field.  The  quantity  Yf 
in  (20)  can  be  measured  using  the  CGS  method,  and  the  question  to  be  considered  is  whether  or 
not  the  right  side  of  this  continued  equality  can  describe  the  measured  angular  variation. 

A  selection  of  CGS  interferograms  obtained  during  dynamic  crack  propagation  in  impact 
loaded,  three-point-bend  specimens  of  PMMA  are  shown  in  Figure  2.  These  patterns  were  recorded 
by  means  of  a  rotating  mirror  high  speed  camera  and  a  pulsed  laser  light  source;  the  details  of  the 
experimental  setup  are  described  by  Krishnaswamy  et  al  (1990).  In  this  case,  the  material  is  trans¬ 
parent  and  the  method  is  used  in  transmission  mode.  The  fringes  in  Figure  2  correspond  to  level 
curves  of  the  quantity  d(cru  + 

In  the  three-point-bend  specimen  of  thickness  h,  the  crack  initiates  abruptly  from  a  blunted 
notch  tup  under  the  action  of  the  stress  wave  loading  provided  by  the  impacting  tip  and,  thereafter, 
it  propagates  at  a  nearly  constant  speed.  Figure  3  shows  the  analysis  of  an  interferogram  at  a  time 
of  t  =  '20  ns  after  crack  growth  initiation.  The  crack  speed  is  approximately  300  m/s.  In  this 
early  time  after  fracture  initiation,  it  is  expected  that  transient  effects  due  to  the  abrupt  initiation 
event  and  associated  stress  wave  propagation  are  still  important.  The  figure  shows  the  variation 
of  Yf  versus  normalized  distance  ri/h  as  given  in  (20)  along  the  radial  lines  corresponding  to 
0t  =  0,30°,  45°,  120°  emanating  from  the  moving  crack  tip.  It  appears  from  the  figure  that  there  is 
no  region  (within  the  framework  of  plane  stress)  over  which  the  measured  Yf  is  constant.  Indeed, 
the  observed  values  scatter  by  as  much  as  several  hundred  percent  about  a  mean  value.  Thus,  there 
appears  to  be  no  region  in  which  a  square  root  singular  stress  intensity  factor  field  dominates  the 
near  tip  field. 

Figure  3  also  includes  the  results  of  applying  a  least  squares  procedure  to  fit  the  expression 
on  the  right  side  of  (20)  to  the  experimentally  measured  values  of  Yf  as  obtained  from  the  inter¬ 
ferogram.  The  procedure  is  based  on  the  values  of  Yf  measured  at  several  points  along  each  radial 
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line.  Because  Yf  and  rj,0;  are  known  at  each  such  point,  (20)  provides  a  linear  equation  for  values 
of  Pk,  k  =  1 ...  ,7  and  Kj  at  each  point.  A  best  fit  solution  of  the  over-determined  system  of  linear 
equations  for  0k  is  then  determined  by  a  standard  least  squares  procedure.  Each  curve  found  in 
this  way  (shown  as  dotted  in  Figure  3)  corresponds  to  the  transient  expansion  on  the  right  side  of 
(20)  and  each  curve  apparently  captures  the  essential  features  of  the  data  for  all  observations  made 
outside  of  the  zone  of  three  dimensional  effects,  namely,  for  r*  >  h/2  (Rosakis  and  Ravi-Chandar 
(1986)).  It  is  emphasized  that  the  curves  all  correspond  to  the  same  values  of  0k,  so  that  the  fit  is 
consistent  in  both  the  radial  and  angular  directions. 

When  the  CGS  method  is  applied  in  reflection  mode,  the  fringes  obtained  correspond  to  level 
curves  of  duz/dxi  where  U3  is  the  out-of-plane  displacement  of  the  surface  of  the  specimen.  For 
deformation  corresponding  to  a  state  of  plane  stress, 

du^  uh  d  (  s 

In  light  of  this  relationship,  the  expression  equivalent  to  (20)  is 


n 


2y/2ZEr*12  du3  _  d  \  mcos|g/  cos  f  0, 

v)vhcos  §0;  cos  ~6t  dxi  1  i[1  cos§0<  2  cos§0/. 


+  ^3(<)rf/a 


cos 


+  r; 


,  cos|0j  ,  .  .  e  .  .cos  lOi1 

1  ^(0 _ J"3~  +  h(t)  +  h(t) 


COS  i  9 1 


COS  \di  \ 


(21) 


where  6k,  k  =  1,...,7  are  functions  of  crack  tip  speed  and  time.  Again,  for  convenience,  the 
left  side  of  (21)  is  described  by  Zf(ri,dj,t).  Just  as  in  the  case  of  Yf  above,  the  quantity  Z*  is 
essentially  constant  in  any  region  in  which  the  fields  are  represented  accurately  by  a  square  root 
singular  stress  intensity  factor  field. 

To  examine  the  role  of  higher  order  terms  in  the  near  field  asymptotic  expansion  in  the  in¬ 
terpretation  of  data  obtained  by  means  of  light  reflected  from  the  sample  surface,  a  second  series 
of  dynamic  fracture  tests  was  conducted.  In  this  series,  the  specimens  were  again  impact  loaded 
three-point-bend  specimens  of  PMMA.  However,  in  this  case,  a  reflective  aluminum  coating  was 
deposited  on  the  front  planar  surface  of  the  specimen  by  vacuum  deposition.  The  experimental 
arrangement  for  reflection  mode  described  by  Krishnaswamy  et  al.  (1990)  was  used  to  observe  the 
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crack  tip  fields.  In  order  to  provide  a  full-field,  visual  demonstration  of  the  agreement  between 
the  six  term  fit  and  the  experimental  data,  a  reconstructed  fringe  pattern  based  on  a  least  squares 
fit  of  the  data  is  shown  in  Figure  4  (discrete  points  on  the  right  side)  superimposed  on  the  actual 
interferogram  (contrast  fringes  on  the  left  side).  In  addition,  the  broken  lines  added  to  the  left  side 
of  the  figure  represent  the  fringe  pattern  reconstructed  from  the  full,  higher  order  least  squares  fit. 
The  higher  order  fringe  pattern  clearly  matches  the  observed  fringes  very  well,  except  close  to  the 
crack  tip  where  the  large  influence  of  three  dimensional  effects  precludes  a  direct  comparison.  For 
purposes  of  comparison,  a  synthetic  fringe  pattern  generated  from  only  the  term  in  the  crack  tip 
field  expansion  associated  with  the  stress  intensity  factor  K1}  is  also  shown  as  the  solid  lines  in  the 
right  half  of  the  figure.  This  pattern  differs  significantly  from  the  actual  pattern,  suggesting  that 
the  higher  order  terms  in  (18)  can  have  a  significant  influence  on  the  crack  tip  stress  field  over  the 
region  in  which  fringes  are  commonly  observed. 

The  discrete  points  in  Figure  5  show  the  variation  of  Zf(rt,0[,t)  with  radial  distance  ri/h 
from  the  crack  tip  for  several  radial  lines  emanating  from  the  crack  tip.  As  in  the  corresponding 
case  of  transmission  mode  observations,  the  quantity  Zj  is  not  constant  over  any  significant  part 
of  the  region  of  observation,  as  would  be  expected  if  the  deformation  field  there  would  be  a  square 
root  singular  stress  intensity  factor  field.  On  the  other  hand,  the  discrete  points  representing 
measured  values  closely  follow  the  curves  obtained  from  a  least  squares  fit  to  the  six  term  asymptotic 
expansion,  at  least  for  rt/h  >  0.5 

To  further  illustrate  the  transient  nature  of  the  crack  tip  fields,  the  time  histories  of  the  stress 
intensity  factor  A'j (<)  and  of  the  coefficients  0k{t)  inferred  from  the  data  during  a  representative 
dynamic  fracture  experiment  are  shown  in  Figure  6.  Each  of  the  quantities  is  normalized  by  its  mean 
value  for  the  entire  test.  Even  though  the  crack  speed  was  constant  during  the  entire  observation, 
the  figure  clearly  shows  the  transient  nature  of  the  crack  growth  process.  For  the  first  50  ns  after 
initiation  of  crack  growth  from  a  stationary  blunted  notch,  there  is  significant  fluctuation  in  the 
values  of  /?*(<),  k  =  2, . . . ,  7.  The  transient  nature  of  a  crack  tip  field  of  this  kind  during  the  early 
stages  of  crack  growth  was  anticipated  in  the  analytical  results  of  Ma  and  Freund  (1986). 
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The  Case  of  Transient  Crack  Growth  at  Nonuniform  Velocity 

The  foregoing  analysis  of  near  tip  fields  during  transient  crack  growth  is  restricted  to  the  class  of 
situations  in  which  the  crack  speed  is  constant.  As  such,  the  analysis  represents  a  full  development 
of  some  ideas  introduced  in  Freund  (1990).  The  case  of  transient  crack  growth  at  nonuniform  speed 
can  be  treated  in  much  the  same  way  but,  as  might  be  expected,  the  details  are  more  complicated. 
In  this  section,  a  three  term  expansion  for  the  first  stress  invariant  is  obtained. 

As  in  the  case  of  constant  speed  crack  growth,  (7a)  and  (7b)  have  solutions 


$o(n,  01,  0  =  Ao(0rf/2  cos  1 9i , 
ri,6l,t)  =  Ai(t)rf  cos 20, . 


(22) 


Up  to  this  point,  the  expressions  appear  to  be  the  same  as  for  constant  speed  growth.  However, 
they  differ  in  the  fundamental  respect  that  the  coordinates  ri,Bi  are  now  functions  of  time.  From 
their  definitions,  it  is  the  crack  speed  that  determines  the  degree  of  distortion  of  these  coordinates 
from  the  physical  polar  coordinates,  and  the  crack  speed  is  now  a  function  of  time. 

For  m  =  2,  the  governing  equation  (8)  takes  the  form 


2  d2$2  d$2 

a‘  drjf  +  di f2 


y/vjt)  d$o 

c(2  drji 


(23) 


If  a  transformation  to  polar  coordinates  (r/,^)  is  made  and  (22)i  is  used,  (23)  becomes 


=  r}/2  cos  \9i  ( Dx  [A0]  +  2 B  sin2  \0,) 


(24) 


where 


Dl[A0\  = 


y/v(t)  d 
a^c2  dt 


4\A(0(l  +  a2)  . 

v^MA(r)  Kl() 


g(|)- 

CfOinAyv)  dt 

The  general  solution  of  (24)  for  mode  I  crack  growth  is 

^2 (r<,0j,t)  =  r5/2  [Aq]  cos  +  B  ^^cos|0(  -  -  cos3  +  A2cos|^ 


(25) 


where  it  is  clear  that  the  last  term  is  the  homogeneous  solution  of  the  equation. 
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A  three  term  asymptotic  expansion  in  powers  of  rt  for  the  first  stress  invariant  can  be  obtained 
by  substituting  the  expressions  for  $o(7/,0/,t),  $1  {rt,0i,t)  and  $2 (r/,0/,t)  into  (10),  with  the  result 
that 

<7ll  +  <722 


3^  cosjg/  2v2 
2 P(cf~c])  4 c2  0  T)n  +  c?  1 


+ 


+ 


f  15V2  \  (  v2  \  v 2 

|^r^2Cos|fl,  +  D'[A0]  |^1  -  cos  \9,  +  cos§0, 

B  [2  (1  _  64^f)  C°S  2°‘  ~  2  i1  ~  2tf)  C°s3  261 

/2  +  o(r2). 


(26) 


—  2cf  C°S  2*'  + 


T“TH}r" 

It  is  readily  verified  that  B  =  0  for  constant  speed  crack  growth,  and  that  this  expression  reduces 
to  the  case  of  (18)  of  equivalent  order.  It  is  noted  that  crack  tip  acceleration  enters  the  expansion 
only  through  the  coefficients  D1[Aq\  and  B. 


Explicit  Expansion  for  the  Broberg  Solution 

In  order  to  illustrate  the  role  of  the  higher  order  terms  in  the  near  tip  expansion  in  a  more  explicit 
way,  the  solution  of  a  particular  boundary  value  problem  concerned  with  elastodynamic  crack 
growth  is  considered.  This  is  the  plane  strain  problem  of  a  crack  growing  symmetrically  from  zero 
initial  length  at  constant  rate  under  uniform  remote  tensile  stress  The  plane  of  deformation 
is  the  x'^Xj-plane  and  the  crack  lies  in  the  interval  —  vt  <  x\  <  1 it,  x'2  =  0,  where  v  is  the  constant 
speed  of  either  crack  tip.  This  is  the  problem  first  analyzed  by  Broberg  (1960). 

An  expression  for  the  first  stress  invariant  directly  ahead  of  the  crack  tips  can  be  obtained 
from  equations  (6.3.43)  and  (6.3.50)  of  Freund  (1990).  On  the  plane  x'2  =  0, 

<Tu  +<722  =  -2<roo/(v/c,)/i  (l-  ^  J  ^-—577^  (27) 

where  h  =  1/v,  I(v/cs)  is  a  known  function  of  v,  and  /(£)  =  (c~2  -  2£2)/(£2  -  c/~2)1/,2(/i  +  £)3/2. 
Focusing  on  the  crack  tip  moving  in  the  positive  indirection,  if  this  expression  is  expanded  in 
powers  of  X!  =  x',  -  vt  near  xi  =  0  then 


<7ll  +  <722  =  W{V) 
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(28) 
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where 


_  2(l  +  a*)(a?-g2) 
[)~  A(«) 


*/(0  = 


c2J(v/cs)&(v) 


v2ai 


3 y/xvt . 


If  the  expansion  (28)  is  then  compared  with  the  general  expansion  (18)  and  terms  of  like 
powers  in  distance  from  the  crack  tip  are  equated,  then  explicit  equations  for  the  coefficients  in  the 
expansion  are  obtained  as 


Ai  =  0 


£>l[/lo] 


15  v 2 

+  ^2  —  — 
4  c\ 


W(v)Kj(t) 

2 Li(cf  jc\  -  \)vt\Z2n 


'1  f'(h)h ' 

.2  +  f(h)  .  • 


(29) 


From  its  definition,  /^[To]  can  be  expressed  in  terms  of  the  stress  intensity  factor.  Consequently, 
the  coefficients  of  the  two  terms  in  (18)  corresponding  to  the  particular  and  homogeneous  solutions 
of  the  differential  equation  for  m  =  3  in  the  asymptotic  expansion  (18)  are 


V‘  -  .  Kj(t) 


(30) 


^2-T  =  92(v) 


where  the  two  functions  of  crack  tip  speed  g\{v)  and  gi(v)  are  shown  in  Figure  7  for  Poisson’s  ratio 
v  =  1/3. 


In  view  of  the  fact  that  the  coefficients  in  (30)  are  proportional  to  t1/2,  it  is  evident  that 
the  third  term  in  the  near  tip  asymptotic  expansion  of  the  first  stress  invariant  can  be  very  large 
during  the  early  stages  of  crack  growth,  possibly  dominating  the  square  root  singular  term  under 
certain  conditions.  This  result  clearly  illustrates  the  strong  bearing  that  transients  may  have  on 
the  question  of  A'^-dominance  of  the  near  tip  field,  even  if  the  crack  tip  speed  is  constant. 
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Concluding  remarks 


Motivated  by  the  wave  propagation  characteristics  of  the  fields  in  the  vicinity  of  a  growing  crack  tip, 
an  approximation  of  the  near  tip  field  in  the  form  of  an  asymptotic  expansion  has  been  introduced. 
The  leading  term  in  the  expansion  of  the  local  stress  field  is  the  familiar  stress  intensity  factor 
contribution,  that  is,  it  is  square  root  singular  in  the  radial  distance  from  the  crack  tip  and  its 
coefficient  is  proportional  to  the  instantaneous  value  of  stress  intensity  factor.  The  higher  order 
terms,  on  the  other  hand,  take  into  account  the  recent  past  history  of  stress  intensity  factor  and 
crack  motion,  in  order  to  reflect  the  transient  nature  of  the  local  field.  These  terms  involve  the 
time  derivatives  of  stress  intensity  factor  Kf(t)  and  crack  speed  v(t).  The  expansion  is  based  on 
an  assumption  that  the  fields  are  indeed  two  dimensional  right  up  to  the  crack  tip. 

For  the  case  of  constant  crack  speed  but  time  dependent  stress  intensity  factor,  the  expansion 
is  given  explicitly  to  six  terms.  The  coefficients  are  evaluated  on  the  basis  of  experiments  involving 
dynamic  crack  growth  in  PMMA  plates  under  the  assumed  conditions.  The  measurements  of  the 
crack  tip  field  were  obtained  by  means  of  the  full  field  Conjugate  Gradient  Sensing  method.  The 
experimental  results,  in  the  form  of  a  gradient  of  the  first  invariant  of  stress  in  the  direction  of  crack 
growth,  appear  to  be  well  described  by  the  expansion.  More  precisely,  over  the  region  near  the 
crack  tip  where  the  elastodynamic  field  is  expected  to  be  essentially  a  two  dimensional  plane  stress 
field,  both  the  observed  radial  and  angular  variation  of  the  data  are  consistent  with  the  higher 
order  expansion. 

The  good  agreement  between  observation  and  the  higher  order  expansion  also  has  a  discon¬ 
certing  aspect.  As  noted  above,  the  experimental  results  and  the  expansion  are  in  close  agreement 
over  a  region  of  significant  size  around  the  crack  tip.  However,  during  the  earlier  phase  of  crack 
propagation,  nowhere  within  this  region  of  agreement  is  the  field  dominated  by  a  stress  intensity 
factor  field.  In  other  words,  even  though  the  specimen  is  large  enough  so  that  reflected  waves  do 
not  influence  the  local  fields,  the  transient  nature  of  the  stress  field  during  the  early  phase  of  crack 
growth  prevents  a  complete  stress  intensity  factor  field  from  becoming  established  outside  the  near 
tip  three  dimensional  zone.  This  feature  erodes  the  value  of  the  stress  intensity  factor  concept 
in  characterizing  the  fracture  resistance  of  the  material.  From  the  discussion  of  the  observations 
in  connection  with  Figures  3  and  5,  it  is  clear  that  a  unique  plane  stress  intensity  factor  can  be 
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identified  in  both  the  data  and  the  model  at  each  instant  of  time.  Furthermore,  Figure  6  shows 
that  the  inferred  value  of  stress  intensity  factor  is  nearly  constant  during  growth,  even  though  no 
fully  developed  stress  intensity  factor  field  is  apparent.  However,  this  stress  intensity  factor  is  not 
strictly  relevant  in  the  spirit  of  fracture  toughness  testing.  The  dilemma  is  particularly  acute  for  the 
case  of  plane  stress  where  a  region  very  close  to  the  crack  edge  is  dominated  by  three  dimensional 
effects,  and  therefore  measurements  in  this  region  cannot  be  used  in  any  two  dimensional  analysis 
of  fields.  The  restriction  is  much  less  severe  for  the  case  of  plane  strain  crack  growth. 

A  re-examination  of  the  experimental  results  of  Ravi-Chandar  and  Knauss  ( 1987)  suggests  that 
this  may  be  an  illustration  of  the  transient  phenomena  discussed  in  this  paper.  In  this  investigation 
the  optical  method  of  caustics  in  transmission  mode  was  used  to  measure  the  stress  intensity  factor 
history  of  a  crack  which  is  loaded  dynamically  by  the  direct  application  of  a  pressure  pulse  on 
its  faces.  In  each  case,  uniform  crack  face  pressure  was  imposed.  Its  magnitude  &o  increased 
linearly  in  time  from  nominally  zero  initial  value  to  a  predetermined  final  value,  with  the  duration 
of  the  increase  always  being  about  T  =  25 /is.  After  initiation,  the  crack  propagated  with  nearly 
constant  speed  v.  A  series  of  experiments  corresponding  to  different  final  crack  face  pressure 
values,  and  therefore  different  rates  of  loading  up  to  fracture  initiation,  were  reported.  Figure  8 
shows  a  comparison  of  stress  intensity  factor  time  histories  obtained  analytically  and  measured 
experimentally  by  caustics.  The  lower  curve  in  Figure  8  demonstrates  good  agreement  between 
experiment  and  theory  both  during  the  dynamic  loading  stage  of  the  stationary  crack  and  during 
crack  growth.  On  the  other  hand,  the  upper  curve  in  Figure  8,  which  corresponds  to  a  higher  rate 
of  loading,  shows  excellent  agreement  of  theory  and  measurements  only  up  to  the  point  of  crack 
initiation.  After  the  crack  started  propagating,  the  experimentally  obtained  stress  intensity  factor 
A'/  overestimates  the  theoretical  result  by  almost  50%.  It  should  be  noted  that  the  delay  time  r 
for  the  crack  to  initiate  after  the  application  of  the  loading  pulse  is  a  factor  of  three  greater  in  the 
experiment  corresponding  to  the  lower  curve  in  Figure  8  than  in  the  experiment  corresponding  to 
the  upper  curve.  The  resulting  crack  tip  speed  in  the  former  figure  is  also  about  half  that  in  the 
latter.  It  is  thus  expected  that  transient  higher  order  effects  of  the  type  demonstrated  above  become 
more  pronounced  in  the  case  of  the  upper  curve.  If  it  is  recalled  that  the  analysis  of  the  optical 
caustics  patterns  is  based  on  the  assumption  that  the  near  crack  tip  field  is  represented  solely  by 


19 


the  r-1^2  singular  term  of  the  asymptotic  expansion  (A'jl-dominance)  then  the  suspicion  arises  that 
caustics  may  not  be  measuring  the  correct  values  of  Kf  after  crack  initiation.  This  becomes  more 
pronounced  in  the  higher  rate  experiments,  where  higher  order  transient  effects  result  in  a  violation 
of  the  basic  assumption  upon  which  the  analysis  of  optical  caustics  is  made. 

An  alternative  illustration  of  the  same  phenomenon  is  seen  in  the  results  of  Krishnaswamy 
and  Rosakis  (1990).  In  this  work  a  “bifocal”  caustics  arrangement  is  used  to  investigate  dynamic 
crack  initiation  and  growth  in  AISI  4340  steel.  At  each  time  instant  during  dynamic  crack  growth 
the  high  speed  camera  used  in  the  experiment  simultaneously  records  two  caustic  patterns  formed 
from  two  initial  curves  at  different  distances  from  the  crack  tip.  Both  caustics  are  then  interpreted 
on  the  basis  of  the  conventional  analysis  with  A'^-dominance.  Since  both  caustics  are  formed  from 
points  surrounding  the  tip  which  He  outside  of  the  near  tip  three-dimensional  zone,  the  values  of 
stress  intensity  factor  obtained  by  bifocal  caustics  pairs  are  expected  to  be  identical  unless  the 
assumption  of  Kf  dominance  is  violated.  The  experimental  results  show  differences  in  Kf  values 
of  up  to  40%.  This  suggests  that  higher  order  transient  terms  of  the  nature  discussed  in  this  paper 
should  be  included  in  the  analysis  of  optical  caustics  in  order  for  the  method  to  furnish  accurate 
values  of  stress  intensity  factors. 

A  long-standing  issue  of  fundamental  importance  in  dynamic  fracture  research  is  the  connection 
between  the  dynamic  fracture  toughness  and  the  crack  tip  velocity.  The  debate,  for  the  most  part, 
has  centered  around  the  question  of  uniqueness  of  a  relationship  between  Kf  and  v.  Kobayashi 
and  Dally  (1980),  Rosakis,  Duffy  and  Freund  (1984),  and  Zehnder  and  Rosakis  (1989),  among 
others,  provided  data  sets  that  seem  to  indicate  that  the  Kf  —  v  relation  is  reasonably  viewed 
as  a  material  property.  On  the  other  hand,  the  results  of  Kobayashi  and  Mall  (1978)  and  Ravi- 
Chandar  and  Knauss  (1984),  based  on  photo-elasticity  and  the  method  of  caustics  in  transmission 
respectively,  suggest  that  there  is  no  such  correspondence.  While  this  difference  may  be  attributable 
to  differences  in  material  characteristics  or  experimental  conditions,  the  fundamental  difficulty  of 
achieving  a  well  developed  stress  intensity  factor  field  in  transient  crack  growth  experiments  may 
also  be  playing  a  significant  role. 

Another  series  of  experiments  leading  to  results  that  have  yet  to  be  explained  are  those  reported 
by  Kalthoff  (1983),  which  seem  to  indicate  that  the  dynamic  fracture  toughness  could  be  specimen 
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dependent.  In  this  case,  as  well,  it  may  be  possible  to  attribute  the  apparent  specimen  dependence 
of  K'j  to  specimen-dependent  differences  in  the  near-tip  region.  Remarks  on  such  observations  are, 
of  course,  speculative  at  this  point,  and  further  study  is  required  to  assess  the  possibilities. 

On  the  basis  of  some  crack  propagation  experiments  in  which  the  optical  method  of  caustics  was 
used,  Takahashi  and  Arakawa  (1978)  proposed  that  the  instantaneous  value  of  dynamic  fracture 
toughness  of  their  material  depended  on  the  instantaneous  crack  tip  acceleration.  As  shown  in 
this  investigation,  however,  the  near-tip  stress  field  expansion  involves  crack  tip  acceleration  in  its 
third  or  higher  order  terms  (see  (26),  for  example).  As  a  result,  caustic  patterns  obtained  from 
regions  where  higher  order  terms  are  important  will  exhibit  acceleration  effects.  However,  if  caustics 
from  such  a  region  are  interpreted  on  the  assumption  of  A*;-dominance  then  it  would  appear  that 
the  instantaneous  value  of  stress  intensity  factor,  and  thus  of  fracture  toughness  of  the  material, 
depends  on  the  instantaneous  acceleration  of  the  crack  tip. 

The  re-examination  of  the  crack  tip  stress  field  for  elastodynamic  crack  growth  under  transient 
conditions  has  led  to  certain  features  of  the  near  tip  field  that  have  the  potential  for  significantly 
altering  the  crack  tip  field  from  the  commonly  assumed  square  root  singular  stress  intensity  factor 
field.  Furthermore,  a  particular  data  set  has  been  examined  in  light  of  the  theoretical  findings, 
with  the  result  that  the  data  are  well  described  by  a  higher  order  representation  of  the  transient 
crack  tip  field.  These  results  offer  hope  of  resolving  some  long-standing  paradoxes  in  the  area  of 
dynamic  fracture,  and  speculation  along  these  lines  has  been  offered  in  this  section.  In  pursuing 
this  issue,  full  field  optical  techniques  such  as  CGS  or  photoelasticity  naturally  are  advantageous 
because  interpretation  of  data  and  thus  measurement  of  Kf  does  not  hinge  on  the  assumption  of 
A'/-dominance.  For  both  techniques  a  higher  order  transient  analysis  of  the  type  presented  above 
can  be  used  to  obtain  the  relevant  coefficients  as  a  function  of  time.  On  the  other  hand,  it  may 
be  possible  to  refine  the  analysis  of  caustics  with  the  hope  that  the  technique  may  exhibit  the 
necessary  sensitivity  to  provide  information  of  the  kind  needed  to  asses  the  importance  of  higher 
order  terms  in  the  crack  tip  field  expansion.  Given  the  extraordinary  experimental  simplicity  of 
the  technique  this  seems  to  be  a  worthwhile  task  for  future  research. 


21 


Acknowledgments 

This  work  was  initiated  during  an  extended  visit  made  by  LBF  to  the  Division  of  Engineering 
and  Applied  Science  at  the  California  Institute  of  Technology.  The  work  at  Brown  University  was 
supported  by  ONR  Contract  N00014-90-J-4051,  and  the  work  at  Caltech  was  supported  by  ONR 
Contract  N00014-90-J-1340. 

References 


Broberg,  K.  B. 

1960 

Arch,  fur  Fisik  18,  159. 

Dally,  J.  W. 

1987 

in  Static  and  Dynamic  Photoelasticity  and 
Caustics  (edited  by  A.  Lagarde),  p.  247, 
Springer- Verlag,  New  York. 

Freund,  L.  B. 

1990 

Dynamic  Fracture  Mechanics,  Cambridge 
University  Press. 

Irwin,  G.  R. 

1960 

In  Structural  Mechanics  (edited 

by  J.  N.  Goodier  and  N.  J.  Hoff),  p.  557. 
Pergamon,  New  York. 

Kalthoff,  J.  F. 

1983 

in  Workshop  on  Dynamic  Fracture  (edited 
by  Knauss,  W.  G.  et  al.),  p.  11.  California 
Institute  of  Technology,  Pasadena. 

Kobayashi,  T.  and 

1978 

Exp.  Mech.  18,  11. 

Mall,  S. 

Kobayashi,  T.  and 

Dally,  J.  W. 

1980 

In  Crack  Arrest  Methodology  (edited 
by  G.  T.  Hahn  et  al),  ASTM  STP  711, 
p.  89.  ASTM,  Philadelphia. 

Krishnaswamy,  S.  and 

1991 

J.  Appl.  Mech.,  to  appear. 

Rosakis,  A.  J. 

Ma,  C.  C.  and 

1986 

J.  Appl.  Mech.  53,  303. 

Freund,  L.  B. 

Ravi-Chandar,  K.  and 

1984 

Int.  J.  Frac.  25,  247. 

Knauss,  W.  G. 

Ravi-Chandar,  K.  and 

1987 

J.  Appl.  Mech.  54,  72. 

Knauss,  W.  G. 

Rosakis,  A.  J.,  Duffy,  J. 

1984 

J.  Mech.  Phys.  Solids  32,  443. 

and  Freund.  L.  B. 

Takahashi,  K.  and 

1987 

Exp.  Mech.  27,  195. 

Arakawa,  K. 

Tippur,  H.  V.,  Krishnaswamy,  S. 

1990 

Int.  J.  Frac.,  to  appear. 

and  Rosakis,  A.  J. 

Yang,  W.  and 

Freund,  L.  B. 

1985 

Int.  J.  Solids  Struc.  21,  977. 

Zehnder,  A.  T.  and 

1990 

Int.  J.  Frac.  43,  271. 

Rosakis,  A.  J. 

22 


Figure  Captions 


Fig.  1  Schematic  diagram  illustrating  the  potentially  complex  dependence  of  the  stress  field  at  a  point 
near  the  crack  tip  on  the  recent  past  history  of  crack  speed  and  stress  intensity  factor  during 
transient  crack  growth. 

Fig.  2  Interferograms  obtained  during  rapid  crack  propagation  in  impact  loaded,  three-point-bend 
specimens  of  PMMA.  The  fringes  are  level  curves  of  the  gradient  of  the  first  stress  invariant 
d(&n  -I-  aii)ldz |,  and  the  times  shown  indicate  the  time  elapsed  since  the  onset  of  crack 
growth. 

Fig.  3  Comparison  of  the  radial  variation  of  the  experimental  measurements  of  Yf(ri,  &i,  t)  as  defined 
by  the  left  side  of  (20)  with  a  fit  to  the  expansion  appearing  as  the  right  side  of  (20)  for  data 
obtained  in  transmission  mode  from  one  PMMA  specimen  (PD- 11).  The  discrete  points  show 
the  results  of  measurements,  and  the  dashed  curves  are  the  result  of  a  least  squares  fit  to  the 
data. 

Fig.  4  Information  on  points  of  constant  stress  gradient  in  the  x\ -direction  obtained  directly  from 
the  photograph  on  the  left  (specimen  PD- 11)  is  shown  in  the  form  of  discrete  points  on  the 
right  side,  and  the  fringe  loops  implied  by  the  points  match  the  actual  loops  on  the  left.  For 
purposes  of  comparison,  synthetic  fringe  patterns  constructed  from  one-term  (solid  lines)  and 
six-term  (dashed  lines)  crack  tip  field  expansions  are  also  shown. 

Fig.  5  Comparison  of  the  radial  variation  of  the  experimental  measurements  of  as  defined 

by  the  left  side  of  (21)  with  a  fit  to  the  expansion  appearing  as  the  right  side  of  (21)  for  data 
obtained  in  reflection  mode  from  one  PMMA  specimen  (PD- 16).  The  discrete  points  show  the 
results  of  measurements,  and  the  dashed  curves  are  the  results  of  a  least  squares  fit  to  the 
data. 

Fig.  6  Time  histories  of  h’j(t)  and  /3fc(t),  each  normalized  by  its  mean  value,  for  the  duration  of  a 
particular  crack  propagation  event. 

Fig.  7  Functions  of  crack  tip  speed  appearing  in  (30)  for  v  —  1/3,  as  determined  by  the  Broberg  crack 
propagation  solution. 

Fig.  8  Dynamic  stress  intensity  factor  versus  time  for  fracture  initiation  and  propagation  in  Homalite- 
100  due  to  sudden  application  of  crack  face  pressure  (Ravi-Chandar  and  Knauss,  1987).  With 
reference  to  the  legends,  op  is  the  level  of  crack  face  pressure  achieved  after  rise  time  T,  and 
the  crack  begins  to  propagate  at  constant  speed  v  after  elapsed  time  r  from  the  instant  that 
loading  begins. 


Fig.  2  Interferograms  obtained  during  rapid  crack  propagation  in  impact  loaded,  three-point-bend 
specimens  of  PMMA.  The  fringes  are  level  curves  of  the  gradient  of  the  first  stress  invariant 
d(au  +  022 ) I dx i ,  and  the  times  shown  indicate  the  time  elapsed  6ince  the  onset  of  crack 
growth. 


Fig.  4  Information  on  points  of  constant  stress  gradient  in  the  zj-direction  obtained  directly  from 
the  photograph  on  the  left  (specimen  PD- 11)  is  shown  in  the  form  of  discrete  points  on  the 
right  side,  and  the  fringe  loops  implied  by  the  points  match  the  actual  loops  on  the  left.  For 
purposes  of  comparison,  synthetic  fringe  patterns  constructed  from  one-term  (solid  lines)  and 
six-term  (dashed  lines)  crack  tip  field  expansions  are  also  shown. 


Normalized  stress  gradient 


Fig.  3  Comparison  of  the  radial  variation  of  the  experimental  measurements  of  Yf(ri,d[,t)  as  dthned 
by  the  left  side  of  (20)  with  a  fit  to  the  expansion  appearing  as  the  right  side  of  (20)  for  data 
obtained  in  transmission  mode  from  one  PMMA  specimen  (PD-11).  The  discrete  points  show 
the  results  of  measurements,  and  the  dashed  curves  are  the  result  of  a  least  squares  fit  to  the 
data. 


Normalized  displacement  gradient 


Fig.  5  Comparison  of  the  radial  variation  of  the  experimental  measurements  of  Zj(r;,  #/,<)  as  defined 
by  the  left  side  of  (21)  with  a  fit  to  the  expansion  appearing  as  the  right  side  of  (21)  for  data 
obtained  in  reflection  mode  from  one  PMMA  specimen  (PD-16).  The  discrete  points  show  the 
results  of  measurements,  and  the  dashed  curves  are  the  results  of  a  least  squares  fit  to  the 
data. 
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Time  histories  of  A'j(t)  and  0k{t),  each  normalized  by  its  mean  value,  for  the  duration  of  a 
particular  crack  propagation  event. 
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Fig.  8  Dynamic  stress  intensity  factor  versus  time  for  fracture  initiation  and  propagation  in  Homalite- 
100  due  to  sudden  application  of  crack  face  pressure  (Ravi-Chandar  and  Knauss,  1987).  The 
lower  and  upper  curves  show  results  for  relatively  low  and  high  final  values  of  crack  face 
pressure,  respectively.  With  reference  to  the  legend,  cro  is  the  level  of  crack  face  pressure 
achieved  after  rise  time  T,  and  the  crack  begins  to  propagate  at  constant  speed  v  at  elapsed 
time  r  after  load  application. 


